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Abstract 

We report a numerical study of the effect of interface alloying and band-alignment on the Auger 
recombination processes of core/shell nanocrystals. Smooth interfaces are found to suppress Auger 
recombination, the strength of the suppression being very sensitive to the core size. The use of 
type-II structures constitutes an additional source of suppression, especially when the shell confines 
electrons rather than holes. We show that "magic" sizes leading to negligible Auger recombina- 
tion [Cragg and Efros, Nano Letters 10 (2010) 313] should be easier to realize experimentally in 
nanocrystals with extended interface alloying and wide band gap. 
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I. INTRODUCTION 



There is current interest in minimizing Auger recombination (AR) processes in semicon- 
ductor nanocrystals (NCs). This is because such processes are considered to be responsible 
for the undesired blinking of NCs, which hinders their use in optical applications.^ When 
an electron- hole pair is generated in a NC by absorption of light, there is a possibility that 
one of the two carriers becomes trapped at the surface. When the next electron-hole pair 
is created, it combines with the remaining carrier. If the remaining carrier is a hole, the 
NC now contains a positive trion (as in Fig. [T^,). If it is an electron, it contains a negative 
trion (as in Fig. |T]d) . At this point there is a competition between radiative electron- hole 
recombination and non-radiative Auger recombination, by which the energy of the recom- 
bining electron- hole pair is transferred to the extra carrier -as shown in Figs. [IJi and [TJd-. 
The excess carrier then moves into a highly excited (typically unbound) state and rapidly 
loses kinetic energy to heat. For optical emission to be free from blinking, the AR process 
must be slower than the radiative recombination one. 

Two promising techniques in slowing down AR processes have been developed in the 
last years: (i) the use of hetero-NCs with radially graded composition,^ and (ii) the use 
of hetero-NCs with quasi-type-II or type-II band alignment Theoretical understanding 
of the AR suppression in (i) was provided by Cragg and Efros™ They showed that the 
smooth confinement potential resulting from the graded composition removes high-frequency 
Fourier components of the electron and hole wave functions, which in turn reduces AR rates. 
Theoretical understanding of the AR suppression in (ii) was provided by us in Ref. [HI where 
we showed that it is originated in the removal of high-frequency components of the carrier 
that penetrates into the shell. 

In this work, we extend Refs. [TOllTTl by performing a systematic study of the influence 
of smooth confinement potentials and band alignment on the AR rates of hetero-NCs. We 
consider not only AR processes involving excess holes (Fig. [T£i), but also excess electrons 
(Fig. [T}d). We find that, as noted in Ref. DUl the softer the confinement potential the slower 
the AR. However, the magnitude of the AR suppression strongly depends on the core sizes, 
with order-of-magnitude differences. A similar behavior is observed when varying the core- 
shell band-offset. In general, moving from type-I to type-II hetero-NC translates into slower 
AR rates, but the actual value is very sensitive to the NC dimensions. The decrease of the 
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AR rate is particularly pronounced when the carrier confined in the shell is the electron 
instead of the hole. Last, we show that the narrow ranges of ("magic") NC sizes leading 
to almost complete suppression of AR reported in Refs. fTUllTol become wider in NCs with 
smooth potential or wide band gap. 




FIG. 1: Schematic of the Auger process under study, assuming the NC contains a positive (a) and 
negative (b) trion. (c) Confinement potential profile used for type-I (left) and type-II NCs (right). 



II. THEORY 

In our calculations, we describe electrons and holes with a one- dimensional, two-band 
Kane Hamiltonian^l 

H= (-& e +y c {x) + E g /2 Kp 

\ R P & h -Vv{x)-E g /2 

where p = —ihd/dx is the momentum operator, m e (h) is the effective mass of the electron 
(hole) disregarding the influence of the valence (conduction) band, hereafter VB (CB). E g 
is the energy gap as defined in Fig. [TJ3, K is the Kane matrix element and V c ^(x) the 
CB (VB) confinement potential. The heigth of V c t v \(x) is given by the core/shell band- 
offset and the shape is depicted in Fig. [IJ:. Note that, at the core-shell interface, V c ( v )(x) 
has a cosine-like profile of width w, which accounts for interfase diffusion. Varying w we 
can monitor the evolution from abrupt interfaces (w — 1 A) to alloyed interfaces spread 
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over several monolayers (w tens of A), accounting for either spontaneous interfase diffusion^ 
or intentionally graded compositionP Equation ([!]) is integrated numerically using a finite 
differences scheme. The resulting electron and hole states have a mixture of CB and VB 
components: 

V't = Xi( x ) Uc{x) + Vi( x ) Uv{x), (2) 

with i = e,h. Here x an d V i u c and u v ) are the envelope (periodic) parts of the Bloch wave 
function corresponding to the CB and VB, respectively. For electrons, the CB component 
X is by far dominant. The opposite holds for holes. 
The AR rate is calculated using Fermi's golden rule: 

2-7T 

= T \ii\V\f)? P{Ef)m- Ef). (3) 

Here, V(r) = l/(e(r + A)) is the Coulomb potential, with e the dielectric constant, 
r = \x\ — x%\ and A a parameter introduced to avoid the Coulomb singularity. Ei and Ef 
are the energies of the initial and final states, respectively. 

To proceed further, we consider the NC are in the strong confinement regime. Let us first 
assume the most relevant AR process is that illustrated in Fig. [l^,. This seems to be the 
case at least in CdSe/ZnS and CdZnSe/ZnSe NCs^H The initial state, |z), is then defined 
by 

l*> = Vfe) ^ (a(l) (3(2) - «(2)) , (4) 

with a standing for spin up (down) projections. The final state, |/), is: 

\f) = ^eVr r h x (x 2 ) + ^r(^i) ^ ni) m - m «(2)) . ^ 

The excited, unbound hole state is taken as ipf^ix) = ^J~j- e lkx u v , i.e. a free electron plane 
wave with negligible component in the CB. L is the length of the computational box. The 
momentum k is determined by energy conservation, k = 2 mo (E® + 2 E® + E g ) /h, where 
mo is the free electron mass, while E® and E® are the single-particle electron and hole energies 
defined with respect to the CB and VB edges, respectively. p(Ef) is the density of states of 
|/), which can be approximated by that of the hole in the continuum (computational box), 
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p(Ef) ~ L/hn ^Jrn^j2~E^ , where E^ x = (k 2 /2m — vbo) is the single-particle energy of the 
plane wave substracting the VB offset. Note that L in p{Ef) cancels the L factor arising 
from the plane wave normalization constants in |(i|l/|/)| 2 . 

Because ipff(x) has no CB component, the Coulomb matrix element in Eq. O) is given 

by: 

(i\V\f) = V2 ((xHVhl V \ Xe e* kX2 ) + (vkVhl V \ Ve e ik ^)) . (6) 

For computational performance and physical insight, it is convenient to rewrite the above 
integrals in the Fourier space. For example, 

1 f°° 

(XhVh\V\x e e ikx2 ) = — J dqV(q)x(q)v(q), (?) 

where V(q) = (e~ tqA E\{— iqA) + e iqA Ei(iqA))/e is the Fourier transform of V(r), with 
Ei(x) = J 00 e~ % jtdt the exponential integral. 14 x(q) = J dx\ Xh{ x i)* Xe(^i) e tqxi , is the 
Fourier transform of the CB components and 77(g) = J dx2T)h{x2)* e^ k ~ q ' X2 that of the VB 
components. 

An analogous development can be carried out for the case of Fig. ^p, where the relevant 
AR process is that involving a negative trion. In such 

(i\V\f) = V2(( VeX e\V\ Vh e lkx >) + ( X eXe\V\ Xh e^)) . (8) 

For the simulations, we take material parameters close to those of typical II- VI NCs: 
m e = 0.15 mo, m h = 0.6 m , E g = 1.5 eV, K = (21eV/2m ) 1 / 2 , e = 6 and A = 10" 2 A. The 
CB (VB) offset at the shell is cbo= 0.5 eV (vbo= 0.5 eV), unless otherwise stated. A high 
potential barrier, 4 eV, is set at the external medium. 

III. RESULTS AND DISCUSSION 
A. Effect of interface alloying 

We start by studying a type-I core/shell NC, where electrons and holes are well confined 
in the core, surrounded by a L s = 20 A shell acting as a barrier. We consider AR involves 
an excess hole (Fig. flk scheme) and model the influence of diffusion around the interface by 
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FIG. 2: (a) Contour plot of the AR rate as a function of core size and interfase thickness in type-I 
NCs. (b) Vertical cross-sections at L c = 25 A (solid line) and L c = 29 A (dotted line). The AR 
process considered is that involving an excess hole. 

shaping the confining potential from abrupt (w = 1 A) to smooth (diffusion spreading over 
w = 30 A) for different core sizes. The resulting AR rates are plotted in Fig. [2^,. As expected 



from simple "volume" scalingp^, in general the bigger the core the slower the AR rate. In 
addition, for a given w, one can see that increasing the core size leads to periodic valleys 
where AR is strongly suppressed (e.g. at L c rj 20, 30, 40 A for w = 1 A). These "magic sizes" 
were first predicted by Efros and co-workers.^ They are originated in destructive quantum 
interferences between the initial and final states of the Auger scattering, and should be 
detectable at low temperature and pumping power We now pay attention to the influence 
of interface smoothness. In general, softening the interface potential reduces the AR by 
up to two orders of magnitude) 10 * 16 ! However this trend is neither universal nor monotonic 
because, as can be seen in Fig. [2^,, softening the potential also changes the position of the 
suppression valleys. This point is exemplified in Fig. [2]d, which illustrates the vertical cross- 
sections highlighted in Fig. [2^, with dashed lines. For some core sizes (L c = 25 A) there is 
a strong decrease of r^ 1 with w. For others, which are initially in a valley of AR (L c = 29 
A), softening the potential provides no benefit up to large values of w (w > 12 A). The shift 
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in the position of the suppression valleys with increasing w is due to the decrease of the 
effective NC size as the shell material diffuses into the core. Interestingly, while the valleys 
are very narrow for w = 1 A (AL C ~ 1 A), they become wider for w = 30 A (AL C ~ 4 A). 
This should facilitate the experimental observation of these minima, since the most precise 
growth techniques rely on monolayer deposition (~ 4 A each monolayer in typical II- VI 
NCs). 
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FIG. 3: Fourier transforms of the CB and VB terms in the Auger matrix element of Eq. (J7]). The 
calculation corresponds to an AR minimum of Fig. [2^,, namely L c = 29 A and w = 1 A. 

Next we provide some insight into the conditions favoring the appearance of magic sizes. 
In the reciprocal space, the minima occur when there is a near cancellation of Eq. [7} Fig. [3] 
shows a typical Fourier spectrum of x(g) and fj(q) at a magic sizeP-^ Note that, while \{q) 
is centered at q = 0, 17(g) is centered at q = k. If the k value is such that one of the periodic 
zeros of 17(g) coincides with the nearby maximum of x(q), as in the figure -see vertical arrow- 
, 17(g) and x(q) are in anti-phase and the AR minimum is obtained. The NC sizes which 
lead to this situation can be approximated with a simple model. Let the potential in Eq. [T] 
be that of a square well with infinite walls and width L c . If we include interband coupling 
by perturbation theory on the basis of the two lowest CB and VB states, the ground state 
wave functions are: 



7T X \ ( 2 7T X 

a, cos I — J u c + b e sin [~^~ 1 "' 

2 7T X \ f 71 X 

a h sin 1 — — \u c + b h cos ( — | //, 



(9) 
(10) 
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where a and b are constants. The Fourier transform of the VB components, fj(q), has zeros 
at qo — k ± j- (2n + 1), with n = 1,2,3. . .. On the other hand, the Fourier transform 
of the CB components, x(q), has the first local maximum at qu L c ~ — 2tt. We can then 
force qM = qo, obtaining k = —tt (2n + 3)/L c . Since k = — a/2 m (E® + 2 + E g ) / ft, using 
particle- in-the-box energies we get: 



m e L\ 



+ 2 



7T^ 



m h L 2 c 

Isolating L c in the above expression leads to: 



+ 



IE, 

h 2 



71 
1~r 



(2n + 3). 



(11) 



tin a/ (2n + 3) 2 — (m /m e + 2m /rrih) 



(12) 



^/2E~m 

For the system in [2^,, Eq. (12) predicts "magic" NC sizes at L = 31, 42 and 53 A, in close 
agreement with the numerical results at w = 1 A. The expression also reveals that wide gap 
materials are more prone to display AR minima. 
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FIG. 4: Same as Fig. J2jja) but for an excess-electron process. 



So far we have considered AR processes involving an excess hole. Very similar results are 
obtained if the process involves an excess electron instead. To illustrate this point, in Fig. [4] 
we reproduce Fig. [2^, but now considering the negative trion process. In general, the AR 
are a few-times slower, but the trends are the same and even the position of the AR minima 
are nearly the same. This implies that magic sizes can suppress AR for either of the two 
process. 
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B. Effect of band alignment 



In this section we study how the core/shell potential height influences the AR rates. 
CB and VB offsets are initially set at 0.5 eV, then we lower the cbo -Fig. [5^i- or the vbo 
-Fig. [5]d-. One can see that lowering the offsets from 0.5 eV to 0.0 eV has a rather weak 
effect on r^ 1 . However, moving to negative offsets (i.e. switching from type-I to type-II 
band- alignment) brings about important modifications. Namely: (i) the position of the 
AR suppression valleys is shifted. This is due to the changes in the effective NC size seen 
by electrons and holes, (ii) The width of the suppression valleys increases, especially for 
negative cbo (Fig. [5^,). (iii) Lowering cbo reduces r^ 1 by orders of magnitude. This is in 
agreement with the experiments of Oron et alP, where using CdTe/CdSe NC -with the 
hole confined in the core and the electron in the shell- lead to a strong decrease of the AR 
rate. Surprisingly, lowering vbo instead (Fig. ^p) barely reduces rj 1 . In other words, using 
CdSe/CdTe NCs would provide much less benefit. 
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FIG. 5: Contour plot of the AR rate as a function of core size and core/shell band-offset, (a) 
variable cbo; (b) variable vbo. The shell thickness is L s = 20 A and w = 1 A. 

The origin of the asymmetric behavior of electrons and holes in Fig. [5] can be ascribed to 
the different effective masses. To clarify this point, in Fig. [6] we compare the density charge 
of the dominant electron and hole wave function components in type-II NCs. Because of 
the lighter effective mass, the electron (xT) displays stronger tunneling across the core. 
The resulting wave function is smoother, which translates into less high-frequency Fourier 
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components than the hole, and hence slower AR. 




FIG. 6: Density charge of the dominant electron and hole wave function components in a type-II 
NC where the carrier is confined in the shell. The gray area represents the confinement barrier. 
L c = 33 A and L s = 20 A. 

We close by noting that in three-dimensional systems, quantitative differences will prob- 
ably arise from our estimates. Yet, the different role of cbo and vbo should persist, as the 
effective masses of the two carriers are still different. Also, we have checked that the weak 
effect of vbo on the AR rate holds for processes involving an excess electron. 



IV. CONCLUSIONS 



Using a two-band Kane Hamiltonian, we have studied how the interface potential shape 
and height influence the AR rates of heteronanocrystals. As previously notecP^, a smooth 
confinement potential -which may result from either spontaneous interfase diffusion or in- 
tentionally graded composition-, reduces the AR rate. However, the magnitude of this 
reduction strongly depends on the core size. This is due to the dependence of the "magic" 
sizes suppressing AR on the interface thickness. Indeed, the range of "magic" sizes increases 
with the interface thickness, which should facilitate their experimental detection. Switching 
from type-I to type-II band-alignment further reduces the AR. For moderate band-offsets 
(fraction of eV), this effect is more pronounced when the shell hosts the electron (instead 
of the hole). This is because the stronger tunneling of electrons enables the formation of 
smooth wave functions. These results are valid for AR processes involving either an excess 
electron or an excess hole. 
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